Surface tension and dynamics of fingering patterns 



F.X. Magdaleno, J. Casademunt 
Departament d'Estructura i Constituents de la Materia 
Universitat de Barcelona, Av. Diagonal, 647, E-08028-Barcelona, Spain 



0\ 



in 

(N 



> 
in 
o 
o 
l> 
o 
l> 

On 
O 

C/3 



Oh- 



We study the minimal class of exact solutions of the 
Saffman- Taylor problem with zero surface tension, which con- 
tains the physical fixed points of the regularized (non-zero 
surface tension) problem. New fixed points are found and the 
basin of attraction of the Saffman- Taylor finger is determined 
within that class. Specific features of the physics of finger 
competition are identified and quantitatively defined, which 
are absent in the zero surface tension case. This has dra- 
matic consequences for the long-time asymptotics, revealing 
a fundamental role of surface tension in the dynamics of the 
problem. A multifinger extension of microscopic solvability 
theory is proposed to elucidate the interplay between finger 
widths, screening and surface tension. 

PACS numbers: 47.54.+r, 47.20.Ma, 47.20.Ky, 47.20. Hw 



The displacement of a viscous fluid by a non-viscous 
one within the gap of a Hele-Shaw cell [jjj, has been one 
of the most studied problems in interfacial pattern for- 
mation for several decades Q . The relative simplicity of 
the problem has made possible an analytical understand- 
ing of the subtle role of surface tension a in the selection 
of the Saffman- Taylor finger |3| , as a prototype of the so- 
called microscopic solvability (MS) scenario of pattern 
selection [Q. More recently, the increasingly interesting 
and controversial issue of the role of surface tension in 
the dynamics of fingering patterns has been raised. 

It is known that the zero surface tension Saffman- 
Taylor (ST) problem is ill-posed as an initial value prob- 
lem and is plagued of finite-time sing ularities (§, §. 
Studies of laplacian growth with zero surface tension, 
however, have proven insightful for instance in cases with 
needle-like growth In the case of smooth interfaces 
which concerns us here, a rich variety of physically rel- 
evant morphologies has been found among solutions of 
the (7 = problem which remain smooth all the time 
(free of finite-time singularities) Q], 0, [||. Given the 
difficulty to get analytical information from the a ^ 
problem, this has raised the question of what part of 
the physics of fingering dynamics, if any, is captured by 
those solutions. Within this spirit, Dai, Kadanoff and 
Zhou explored via numerical simulation the qualitative 
differences of distinct classes of initial conditions || . For 
the so-called pole-like class and for finite-time, the prob- 
lem was concluded to be qualitatively similar with and 
without surface tension. More recently, Siegel and Tan- 
veer [fiof have shown that the regularized problem (i.e. 
vanishingly small a) may differ significantly from the ide- 



alized problem (a — 0) in order one time, and therefore, 
smooth time evolving solutions of the a = problem do 
not coincide, in general, with the limit of solutions of the 
regularized problem. Evidence for this is shown in the 
particular class of single-finger configurations. This re- 
sult, however, does not preclude other situations where 
this might no be the case More generally one can 

find situations where the evolution in the two cases is 
qualitatively equivalent in the sense that a (small) quan- 
titative difference between the two remains bounded for 
all time [Q. Further physical insight is thus still neces- 
sary to clarify the phenomenology which may be appro- 
priately captured by the idealized problem, particularly 
concerning the long-time asymptotics. 

From a different perpective, the MS scenario itself has 
been questioned recently by results of Ref. |f3| where it 
is claimed that in a wide class of nonsingular exact solu- 
tions, the a — dynamics leads naturally to the solution 
predicted by selection theory, without invoking surface 
tension to explain selection. This seems to support the 
claim that surface tension is indeed unessential to the dy- 
namics. In this letter we sustain the opposite conclusion 

& 

Our approach here consists of identifying specific dy- 
namical features which can be viewed as essential to the 
process of finger competition from a physical standpoint, 
and check them in exactly solvable zero surface tension 
cases. This will allow for a precise diagnosis on the phys- 
ical content of the idealized a = problem. 

The evolution equation for the time-dependent confor- 
mal mapping f(w,t) of the interior of the unit circle in 
the complex plane w into the region occupied by the vis- 
cous fluid in the physical plane z = x + iy, in the case 
of zero surface tension and in the co-moving frame (the 
frame moving with the mean interface velocity) can be 
written as 



Re(wd w f{w,t)[l + dtf*{w,t)]) = -1. 



(1) 



A channel of width 2?t is considered along the x axis, and 
periodic boundary conditions are assumed in the y axis 
for simplicity. 

The basic idea is to find a simple, low-dimensional non- 
singular class of solutions of Eq.(Q) which contains the 
physical fixed points of the regularized problem, and com- 
pare the phase space flow topology in both cases. The 
key point is that we do not need to know the exact phase 
space trajectories of the regularized problem, which are 
particularly difficult to obtain even numerically for long 
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times, but only the phase space flow topology. The lat- 
ter can indeed be inferred unambiguously from existing 
empirical evidence both experimental and simulational. 

A class of solutions which satisfies the above require- 
ments is of the form 



f(w,t) = - lnw + d(t) 
+ (1 - A) (ln(l - a(t)w) + ln(l + a*(t)w)) 



(2) 



where a(t) — a! if) + ia"(t). This corresponds generi- 
cally to two unequal fingers. The interface shape has two 
symmetry axes along the fingers, separated a distance w. 
Such symmetry simplifies the analysis but does not affect 
the competition of fingers in any fundamental way fl6fl . 

The case a'{t) = corresponds to a single finger with 
the asymptotic ST shape for a" — > 1. The case a"(t) — 

corresponds to two identical fingers which tend to a 
doubly degenerate ST solution as a' — » 1. For \a(t)\ -C 

1 the ansatz (|^) describes sinusoidal perturbations of a 
planar interface Jl7| . A is a constant of motion and takes 
real values in the interval [0, 1] . For well developed fingers 
A is the total filling fraction of the channel occupied by 
the invading fingers. 

For a = this ansatz is exactly solvable in the sense 
that if we insert Eq. (||) into Eq. ([!]) we obtain a closed set 
of ordinary differential equations for the parameters a', 
a" and d. According to Ref. this case is free of finite- 
time singularities. The parameter d accounts for a global 
displacement and is irrelevant for the present discussion. 

For the sake of discussion and visualization, we find 
convenient to parametrize the phase space in terms of the 



variables u = 1 — a" 2 and r 
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Thus the phase space is the cube [0, 1] x [0, 1] x [0, 1] in 
the (it, r, A) space. In these variables, the time evolution 
is given by the equations 

/ . 3r — 4 — qr(l — ru) 

u = 2ru(l-« — ' 3 

l+gT g (u,r) 

j. 2r(1 r f r 2 (! + ™) + .9(1 - ~ r) (4) 



1 + gT g (u,r) 



where 



T g (u, r) = (1 - g)(2r + g(2r - 1)) - -(1 - gfru 

-gur 2 (l+g(ru-3)) (5) 

and where g = 1 — 2A = const. 

In order to compare with the physical case of a ^ 
we introduce the following construction. Consider a one 
dimensional set of initial conditions (t = 0) of the form 
Eq.(||) surrounding the planar interface (PI) fixed point 
u = 1, r = 1, for a fixed A. We take them infinitessimally 
close to PI in such a way that the interface is in the linear 
regime fl7| . The time evolution from t = — ooto£ = oo 
of this set spans a compact two dimensional phase space 



(it', r') imbeded in the infinite dimensional space of inter- 
face configurations. For a finite surface tension, we know 
from all existing empirical evidence that the above sub- 
space must contain three fixed points, namely the (un- 
stable) planar interface (PI), the (stable) ST single finger 
(1ST) and a saddle fixed point corresponding to the de- 
generate double ST finger (2ST). For finite a (and in the 
high viscosity contrast limit |16|), the 1ST fixed point is 
known to be the universal attractor of this problem so all 
trajectories start at PI and end up at 1ST. The 2ST fixed 
point has a lower dimension attracting manifold defined 
by a" = and will govern the dynamics of finger compe- 
tition. We may define the space (u',r')o as the limiting 
case of a — > (taken after the limits t — > ±oo). 

Since in the linear regime the regularized problem for 
vanishingly small a converges regularly to the a = solu- 
tion, the manifolds (u, r) and (it' , r') must be tangent at 
the PI fixed point, u = 1, r = 1 (see Fig.l). Furthermore, 
(u, r) and (u',r')o must intersect at 1ST and 2ST which 
occur respectively at u = 0, r = 1 and u = l,r = 0, as 
seen directly from Eqs.(|)-(|). Although (0, 1) and (1, 0) 
are fixed points for any A, according to selection theory 
the intersection will occur at A = 1/2. 

Following the topological approach of Ref. |Q, it is 
useful to consider the stream function if), defined as 
the imaginary part of the complex potential <f>(w,t) — 
—f(w, t) — In w, in the co-moving frame. Along the inter- 
face, tp is then a periodic function which provides a nat- 
ural definition of individual aerial growth rate of fingers, 
which we call Atpi, and Alps for the longer and shorter 
fingers respectively. In this simple case these are given as 
maximum-to- minimum differences of the stream function 
extrema along the interface, ip may have only one maxi- 
mum even for two-finger configurations, in which case we 
take Alps = & n d qualify the finger as 'non-growing'. In 
our case, and for finite width fingers, Atp > (growing) 
and = (non-growing) correspond respectively to 
positive and negative tip velocities relative to the mean 
interface position. 

The physical scenario of finger competition which we 
want to test, extracted from experiments and simulations 
|l6| can be briefly described as follows. In the linear and 
early nonlinear regimes two different fingers grow with 
both AipL(t) and Aips(t) increasing with time. When 
the fingers are well developed and the 'growth' function 
G(t) = A *M«)+Afe(t) (with 1-A = AtP l (oo) + AtPs(oo)) 
is of order one, the competition takes over. This is sig- 
naled by an enhanced growth of the 'competition' func- 
tion C(t) = A-fc(tj+Afe(tj as starts to increase at 
the expense of Alps. Existence of competition can thus 
be identified with Alps decreasing with time. The compe- 
tition may be termed 'succesful' when Alps — > asymp- 
totically (C(oo) = 1), that is when a 'growing' finger is 
turned into 'non-growing' due to the presence of another 
finger. As discussed in Refs. pq|, this dynamical elimina- 
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tion of the small finger is associated to topology changes 
in the physical velocity field, which occur via the crossing 
of topological defects through the interface. 

We now discuss the zero surface tension dynamics of 
our ansatz (^), by analyzing the phase portrait of the 
dynamical system denned by Eqs.(||)-(||). This is plotted 
in Fig.l for A = 1/2. In this case Eqs.(||)-(||) can be 
integrated analytically. Dynamical trajectories are of the 
form 



2„,2 



2u — 3ru + r u 
^/u(l — r)(l — ru) 



const 



(6) 



and are plotted as solid lines with arrows. The dashed 
lines are A-independent boundaries. The short-dashed 
line separates the one-finger (above) and the two-finger 
(below) regions. The long-dashed one is the defect 
boundary separating the no-defect (above) and the 1- 
defect (below) regions, with Alps — and Alps > 
respectively. 

Our central result is that the topological structure of 
the phase portrait for a = (Fig.l) is radically differ- 
ent from that of finite surface tension. Its most salient 
feature is the fact that the ST single-finger solution has 
a limited basin of attraction. Part of the flow evolves 
towards a continuum of (attractive) fixed points r = 0. 
The separatrix of the two regions is a critical trajectory 
ending in a new (saddle) fixed point located at u* = 
and r* = 2A/(1 + A). 

For arbitrary A, the line r = is a continuum of sta- 
tionary solutions with coexisting unequal fingers (differ- 
ent widths Ai, A2 with Ai + A2 = A) advancing with the 
same velocity and with tip positions separated in the re- 
direction by a finite distance A = (1 — A) log ^j^=I . 
Solutions of this type have been reported previously 0, 
H . We would like to call the attention upon the fact that 
the screening of the laplacian field, as the mechanism usu- 
ally invoked to explain competition, applies actually to 
the aerial growth of the small finger Alps, which is re- 
duced indeed by the presence of the longer neighboring 
finger, but not to the velocity, which may be in fact the 
same. Concepts such as screening length and time Q 
are therefore meaningful only if the widths of competing 
fingers are constrained to be equal. 

The fixed point (u* , r* ) corresponds to a new type of 
asymptotic stationary solution of the a = ST prob- 
lem. It consists of two fingers with unequal positive 
velocities. The length ratio of the fingers satisfies Jlf| 
limi_oo Ls/Ll = 1/3 independently of A. For the r = 
solutions we have Hindoo Ls/Ll = 1 while, for the 1ST 
fixed point when approached from the two-finger region 
we have lim^oo Ls/Ll = 0. In the latter case, the resid- 
ual non-growing finger which subsists is reminiscent of 
the 'frozen' fingers observed in real experiments. 

Our central point is that, according to the above dis- 
cussion, the possibility of successful competition is asso- 



ciated to the fact that dynamical trajectories cross the 
defect boundary from below (annihilation of topological 
defects [^6|). In Fig.l we see that, for A = 1/2, there 
is no successful competition whatsoever since the critical 
trajectory is located above the defect boundary pl|. 

The cases of A ^ 1/2 and a = are not directly rel- 
evant to the viscous fingering problem, but may be rel- 
evant to other generic situations of laplacian growth in 
the spirit of Refs. ||, and will be discussed in detail else- 
where jl5| ]. Here we will just remark that r*(X) is mono- 
tonically increasing, between r*(0) = and r*(l) = 1. 
Therefore, the basin of attraction of the single finger so- 
lution is larger for narrower fingers. Furthermore, since 
the defect boundary is independent of A, there exists a 
critical A c = 1/3 for which r* crosses the defect bound- 
ary. This implies that, for A < 1/3, there are dynamical 
trajectories which cross the defect boundary from below, 
and therefore the competition is then successful for some 
finite region of phase space. 

In summary, from the analysis of the a — dynam- 
ics of the class Eq. (||) we conclude that (i) only a small 
A-dependent part of phase space behaves qualitatively 
as the a ^ problem, leading to a ST single finger 
(with maybe a residual non-growing finger); (ii) dynam- 
ical elimination of growing fingers does not occur for fin- 
ger widths relevant to the problem of viscous fingering 
(A = 1/2); (iii) the picture of competition based solely on 
laplacian screening is insufficient, since relative widths of 
fingers and not only relative tip positions come into play. 

In order for the 'screening' picture to be valid, an ad- 
ditional dynamical constraint is required to force the fin- 
ger widths to be equal. In growth processes based on 
aggregation of particles, the finger width may be fixed 
by particle size (set to zero in Refs. ||). In the problem 
of viscous fingering such constraint is supplied precisely 
by surface tension. This suggests that an extension of 
MS, which is essentially a static theory, to multifinger 
configurations, may shed new light on the dynamics of 
the problem. 

The generalized multifinger MS scenario can be 
sketched as follows. For two-finger configurations, there 
exists a two parameter continuum family of steady state 
solutions which we can parametrize by A = Ai + A2 and 
p = Ai/(Ai + A2). The cases p = 0, 1 correspond trivially 
to the single-finger case. The case p = 1/2 (two identical 
fingers) is also reducible to single-finger MS in a chan- 
nel of half width. Most interestingly, one can show H 
that, for nonzero surface tension, nontrivial stationary 
solutions with unequal fingers (p 7^ 1/2) exist. In this 
case, surface tension selects an infinite set of values of A 
which differs from the single-finger case, but which scale 
also as (A - 1/2) ~ cr 2 / 3 . Further more, for any given A 
of the above discrete set, there exists another countably 
infinite set of possible values of p with (p — 1 /2) ~ icr 1 / 3 
p9[ . This new set of fixed points of the problem with sur- 
face tension and its physical relevance will be discussed 
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in detail elsewhere |fj]]. It is reasonable to expect that 
all two-finger solutions except the p = 1/2, A = A m (cr) 
(with A m the minimum value of A from MS) lie outside 
the space (u',r'). Despite the fact of being globally un- 
stable, the p — 1/2, A = A m fixed point is then the phys- 
ically relevant one to describe finger competition since it 
has an attracting manifold which includes the PI fixed 
point (linear regime). The process of finger competition 
can thus be pictured as follows. From the linear instabil- 
ity a given number of fingers emerge. As far as this early 
stage is dominated by the most linearly unstable mode (in 
the limit of weak white noise on PI, the emerging con- 
figuration will indeed be nearly periodic), the interface 
is relatively close to the attracting manifold of the nST 
fixed point (n equal fingers with A, = \ m (a)/n). The 
fingers tend thus to adopt the same fingertip curvature 
and select their widths at early stages of the nonlinear 
regime according to single-finger MS theory. The nST 
fixed will then govern the process of competition in the 
sense that the path connecting typical initial configura- 
tions with the single finger attractor, must necessarily 
pass near that saddle point. The phenomenon of com- 
petition is then viewed as the crossover to the unstable 
directions of the nST, A = A m fixed point. Such crossover 
is what is missed in the a = problem, since the unsta- 
ble direction of the 2ST fixed point becomes infinitelly 
marginal (a line of fixed points) in that limit. In the 
terminology of dynamical systems, this reflects the fact 
that Eqs.(||)-(||) are structurally unstable ^0|. We thus 
conclude that surface tension plays a fundamental role 
in the dynamics of finger competition and that, for the 
long time asymptotics, it can only be treated as a 'regu- 
lar' perturbation in a very limited region of phase space 
which excludes multifinger configurations. 
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Fig f . Phase space flow of the dynamical system de- 
fined by Eqs.(|])-(||) for A = f/2. See explanation in text. 
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